One dimensional drift-diffusion between two absorbing boundaries: application to 

granular segregation 
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Motivated by a novel method for granular segregation, we analyze the one dimensional drift- 
diffusion between two absorbing boundaries. The time evolution of the probability distribution and 
the rate of absorption are given by explicit formulae, the splitting probability and the mean first 
passage time are also calculated. Applying the results we find optimal parameters for segregating 
binary granular mixtures. 



I. INTRODUCTION 



The diffusion phenomena had been one of the most 
intensively studied fields in statistical physics. A num- 
ber of textbooks have been written in this field, and in 
many of them the one dimensional case is discussed in de- 
tail . Still, according to our knowledge, the special 
case when the diffusing particle is between two absorbing 
boundaries, and a constant external force is applied, has 
not been completely solved. 

The need for the solution of this problem arose when 
we investigated the motion of a granular particle in a 
vertically vibrated ratchet by means of a 2D computer 
simulation (see Fig. |l]). 




FIG. 1. The setup of the one particle simulation. The two 
dimensional box has an asymmetrical sawtooth shaped base, 
which is sinusoidally vibrated in vertical direction with am- 
plitude A and frequency /. The shape of the sawtooth can 
be described by three parameters: width w, height h, and 
asymmetry parameter a, which is the ratio of the projection 
of the left edge onto the base and w. The particle has three 
parameters: radius r, coefficient of restitution e and friction 
coefficient /j,. The boundary condition can be either periodic, 
reflective, or absorbing. 

For a detailed description of the setup see ||. We 
found that the horizontal motion of one particle can be 
well approximated as drift-diffusion (see Figures |] and ||) , 
and the parameters of the diffusion process, i.e. the av- 
erage velocity and the diffusion constant, depend on the 
parameters of the particle and the ratchet |(| . According 
to our results, it is possible that the average velocities 
of two kinds of particles in the same system have oppo- 



site directions. If the boundaries are open (i.e. absorbing 
boundary condition is applied) , then this setup is capable 
of segregating a binary granular mixture of these parti- 
cles provided that the load rate is chosen so that only a 
few particles are in the system at one time (in this case 
the interaction between the particles can be neglected, 
and we do not have the problem that the transport ve- 
locity also depends on the number of particles) . We need 
the theoretical description to predict e.g. the quality of 
the segregation or the highest possible load rate, and to 
further improve the quality by choosing optimal parame- 
ters, such as where to load the granular mixture in along 
the box. The results of the granular binary mixture seg- 
regation will be published elsewhere. 
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FIG. 2. Simulation result: the probability distribution of 
the horizontal position of a particle which is started from 
x = at zero time, in an infinitely wide system (natural 
boundary condition). The dashed curves show the theoreti- 

cal prediction -j==e tot , which is the probability dis- 
tribution of a drift-diffusing particle in one dimension. The 
velocity v and diffusion constant D are determined by line 
fitting: (x(i)} = vt and ait) 2 — 2Dt (see inset), where 
a(t) 2 = (x(t) 2 ) - (x(t)) 2 . The fitted values are v = 1.42 

2 

and D — 4.42 The parameters of the simulation are the 
following: A — 2 mm, / = 28 Hz, w = 6 mm, h = 8.5 mm, 
a = 0.07, r = 1.2 mm, e = 0.45, and fj, = 0.1. 
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FIG. 3. The same as in Fig. B except for that the resti- 
tution coefficient of the particle is e = 0.6. The fitted diffu- 
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sion parameters are v = —1.60 ^jp and D = 11.43 Note 
that the direction of the velocity is the opposite, which makes 
possible the segregation of a granular mixture consisting of 
particles with e = 0.45 and e = 0.6. 



II. THE DIFFUSION EQUATION WITH BIAS 

The dynamics of the diffusing particle is characterized 
by the diffusion constant D and mean velocity s / 0, 
which is the result of an external field. (The results for 
the v = case will be presented later.) The diffusion 
equation then reads dtp(x,t) — Dd^.p(x ) t) — vd x p(x,t). 
The relevant quantities are D and v, which define a char- 
acteristic time t* = D/v 2 and a characteristic length 
I* = D/v. We use t* and /* to non-dimensionalize the 
problem, all quantities are dimensionless in the rest of 
this section. The dimensionless diffusion equation is 



d t p{x, t) = &lp(x, t) ~ d x p(x, t) 



(1) 



with initial condition p(x, 0) = Po(x) for < x < L, and 
boundary condition p(Q, t) — p(L, t) = for t > 0, where 
L is the system width, and we start at t — 0. The absorb- 
ing boundary condition is equivalent to that the proba- 
bility is zero at the boundaries || . Let us define operator 
F = dl - d x with V F = {(j>(x)\(f>(0) = <j)(L) = 0} n V 2 , 
where T> 2 is the set of twice differentiable functions. If 
4>\ G T>p is an eigenfunction of F with eigenvalue A, 
then exp(Xt)4>\ is a solution of (|l|). If we had an or- 
thonormal and complete eigensystem of F, we could give 
the solution of ([[]) immediately. It is a problem that 
F is not Hermitian, as far as the usual scalar product 
Jq 4>i(x)(j)2(x)dx is concerned. However, F is Hermi- 
tian if the scalar product is calculated with kernel func- 
tion e~ x : (0i 1^2) = /(f e~ x (j)\(x)(j)2{x)dx. It is useful 
to write F in the form F = e x d x e~ x d x . The verifica- 
tion of that F is Hermitian with scalar product (-|-), i.e. 
{4>x\F4>2} = (F(f>i\4>2) is straightforward: the effect of 
operator F can be shifted off from cf>2 to <pi using two 
consecutive partial integration, during which the surface 



terms disappear, because 4>i(x) and fofa) are zero on 
the borders. Since F is Hermitian with scalar product 
(■|-), it has a complete orthonormal eigenfunction sys- 
tem with real eigenvalues. It is straightforward to ver- 
ify that the eigenfuctions are Pk(x) — \J^ eX ^ 2 sin (j^x) 

with eigenvalues = - (^) 2 - \ (fee N+). The solu- 
tion of (Q), if p°(x) is the initial probability distribution, 

00 

is p(x,t) — (p°\Pk)Pk(x)e Xkt . In the rest of the paper 

k=l 

we investigate the case when the initial distribution is a 
Dirac-delta function at aL (0 < a < 1). The weights 

are then (5(x — aL)\pk) = \J~^; e ~ £ ^~ sin(£;7ra), and the 
probability distribution is given by 



p(x,t) 



2^ 

1 x — aL 

= 2L 6 



sin(fc 



k=l 



kit 



X \ 7T . . 

a -I )r z{t) 



-[(¥) 2 +l> 



X \ 7T 

Z 2'^ 



(2) 



where we introduce the notation z(t) = e and the 

00 2 

theta function ^(r, q) = 1 + 2 ^ cos(2rk)q k to obtain 

a closed form J7J. An important quantity is the rate of 
absorption, i.e. the probability current at the borders. 
Using the notation $' 3 (r, q) = d?i9 3 (f, q)\r =r , the proba- 
bility current j{x,i) = —d x p(x,t) + p(x,t) at the right 
and left end is 



^(t)=j(L,t) = ^e— *— 



j^(t) = -j(0,t) 



2ZX 



3 7?'o 



ir[a + 1] 



-,z(t) 



z(t) 



(3) 



With the minus sign in its definition, j<—(t) > 0. 




FIG. 4. The time evolution of the probability distribution, 
when L = 2 and a = 0.3. Note the difference between the 
probability distribution here and in Figs. |^ and |^, which is 
due to the different boundary condition. Inset: the probabil- 
ity current at the left and right boundary. 
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The splitting probablility, i.e. the probability that the 
particle is absorbed finally by the left or the right bound- 
ary: = L j^(t)dt, and n_» = L j^(t)dt, and ob- 
viously n<- + = 1 . Using the integral formula 



7T 

2' 6 



ds = 



2l 2 smh(xy/a) 
7r sinh(^Y / o) 



which holds for \x\ < I |g, we get 



-oi 



-ai 



1-e" 



and n_ 



1-e" 



(4) 



(5) 



Another important quantity is the mean first passage 
time, i.e. the average time it takes for the particle to be 
absorbed by any of the boundaries. The probability dis- 
tribution of this time is just the total probability current 
at the boundaries: j ou % = j<_ + The mean first pas- 
sage time is then r = J °° tj ou t(t)dt. The calculation of 
r is straightforward using the derivative of (^) w.r.t. a, 
and the result is 



L(n_> — a). 



(6) 



There is another way to calculate n_> and r, in which 
there is no need for an explicit formula for the time de- 
pendent probability distribution or the current at the 
borders An ordinary differential equation can be 

written for n.— which, using our notation, is —Ld a n^ + 
d^ri-t = with boundary condition n^| Q= o = and 
ri— i\a=i — 1- By direct substitution one can verify that 
(||) is the solution. A similar equation can be written for 
the mean first passage time: —Ld a r + d^r = —L 2 with 
boundary condition r\ a= Q = and r| a= i = 0. It is easy 
to check that m) is the solution for this equation. 



III. DISCUSSION OF THE RESULTS 

First we summarize our results with dimensionalized 
parameters. The probability density is 



p(x,t) 



V v(2l a L-a:] + vt) 

-e -id J 



2DL 

-03 



a 



L 



a - 



(7) 



case (0) and (||) are valid with v — substitution, while 
for the splitting probability and the mean first passage 
time we can take the v — > limit and obtain = a 
and r = a(l- a)L 2 /(2D). 

Since the splitting probability is the most important 
quantity as far as the segregation is concerned, we ana- 
lyze it in detail. Its dependence on all four parameters v, 
D, a, and L can be seen in Fig. ||. It can be easily under- 
stood that with larger (positive) velocity or smaller dif- 
fusion constant the right arriving probability is closer to 
1. It is also trivial that if the starting position is closer to 
the right end, then is larger. However, to understand 
the dependence on the system width, we have to recall 
that in the case of drift-diffusion with natural bound- 
ary condition (i.e. when the system is infinitely wide) 
the expectation value of the position is proportional to 
the time, while the dispersion is proportional only to the 
square root of time. Although this is not exactly true in 
the case of absorbing boundaries, the tendencies remain 
the same, therefore, in case of v > 0, the larger the sys- 
tem width is, the less part of the probability is absorbed 
by the left boundary, and as a consequence, the rest is 
absorbed by the right boundary. 




-1 -0.5 0.5 1 1.5 2 2.5 3 
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FIG. 5. Dependence of the splitting probability n_> on all 
the parameters. The default values for the parameters are 
v = 1.42 £S, D = 4.42 2ai (the fitted values in Fig. §), 
a — 0.5, and L = 20 cm. 



IV. AN APPLICATION: SEGREGATION 



with the notation z(t) = e n Dt / L . The currents at the 
right and left borders are 



. . , 7TV t,(2[o-i]i,+ t ,t) f\a— l]ir , . 
= ^7I e ^ *s ( l —^ J -,»(t) ) - 



2DL 2G -^MT'*<*> 



(8) 



I _ e -aLv/D 

The splitting probability is = — > an< ^ the 

mean first passage time is r = L(n_> — a) /v. In the v = 



As mentioned in the Introduction, our motivation for 
the above calculations has been the application for segre- 
gating binary granular mixtures. For this purpose, now 
let us investigate the segregation properties of the sys- 
tem. 

We begin with an illustrative example, the case of sym- 
metric initial condition, a = 1/2. Denoting v/D by u, 
the splitting probabilities are 



1 



1 



1 + e uL / 2 ' 



1 + e - uL l 2 ' 



(9) 
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For definiteness, let us have a particle with u > 0, for 
example. For a small L, both probabilities are around 
1/2: 



1 



In A 



1 



1 



l/2-u£/8, S3 l/2 + uL/8. 



For large is, 



!_ e -«i/2 



(10) 



(11) 



they tend to zero and one, respectively, in an exponential 
way in L. We can see that, for a small L, it is the diffu- 
sion, the left-right symmetric effect, that determines the 
left and right probabilities. On the other side, for large 
Ls, the drift becomes the dominating effect, it drives the 
particle to the direction of u with a probability that dif- 
fers from 1 only by an exponentially small amount. This 
latter phenomenon makes it possible to use the system 
for segregation. 

To study the segregation properties in more detail, now 
let us have two types of particles, with u\ < < u 2 - We 
neglect the interaction between particles, which is a good 
approximation at low particle numbers. The expectation 
is that the particles with u\ will tend to move to the 
left end and the others to the right end. Therefore, we 
can characterize the quality of the segregation with the 
quantity 



(12) 



Let us put the question whether, for a given L, it is pos- 
sible to choose an optimal a, i.e.,where q is maximal [[)]]. 

Investigation of the function q(a) shows that it indeed 
has one maximum, at 



In 



"l(e" 



-1) 



^optimal 



= l - 



«2<e" 



-1) 



(u 2 - ui)L 



(13) 



For an illustration of the typical dependence of q on a, 
see Fig. |(| 




0.25 
FIG. 6. Dependence of q on a, for u\L 



For practical purposes, we are interested in large Ls. 
In this asymptotic region, with the notation A = U2/\ui\, 



^optimal 



1 + A 1 + A 1 + A 



O 



(14) 



Consequently, 



— 

n^(ui) ss 1 — c\e !+ A 



\ Ul L\ 



n^(u2) S3 1 — c 2 e i+ 



(15) 



, —A 

where C\ = X 1 + x and c 2 = X 1 + x . Therefore, we can 
see that, for large Ls, it is possible to choose such an 
optimal a that both segregation probabilities are expo- 
nentially close to 1 (as functions of L), the quality of the 
segregation is exponentially close to the perfect. 



V. SUMMARY 

We presented a complete solution to the problem of one 
dimensional drift-diffusion (with constant external force) 
between two absorbing boundaries: the probability dis- 
tribution, the probability current at the boundaries (i.e. 
the rate of absorption), the splitting probability and the 
mean first passage time were calculated. The results were 
applied to predict the quality of granular segregation in 
a vertically vibrated ratchet. We found that if the com- 
ponents have opposite drift velocities, the quality of the 
segregation of a binary mixture increases rapidly with 
increasing system width, and as a limiting case, perfect 
segregation can be achieved. Furthermore, when the sys- 
tem width is fixed, we found the place where the granular 
mixture should be loaded into the system to obtain the 
best segregation quality. 
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